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Abstract 

In the WKB approximation the V 2 5 term in Schrodinger's equation 
is subordinate to the (VS 1 ) 2 term. Here we study an anti-WKB ap- 
proximation in which the V 2 S term dominates (after a guess for S 
is supplied). Our approximation produces only the nodeless ground 
state wavefunction, but it can be used in potential problems where 
the potential is not symmetric, and in problems where there are many 
degrees of freedom. As a test, we apply the method to potential prob- 
lems, including the hydrogen and helium atoms, and to </> 4 field theory. 
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1 Introduction 



In quantum mechanics the wavefunction is sometimes written in the form 

tp(r) = exp S(r). (1.1) 
Schrodinger's equation then takes the form 

V 2 S + \VS\ 2 = -p- [V(r) — E\ . (1.2) 

The WKB approximation ensues when the term V 2 S* is dropped in leading 
order and later incorporated as a correction. It is commonly understood 
that this step is justified in the semi-classical regime where the gradient of 
the de Broglie wavelength, VA, has a magnitude much smaller than one, and 
the wave function oscillates many times over distances that characterize the 
variation of the potential. 

In this paper we consider the alternate ordering of terms where the term 
\VS\ 2 is dropped in leading order and later treated as a correction. This is 
an anti-WKB approximation. As we shall see, the anti-WKB approximation 
is valid deep in the quantum regime where the wavefunction is nodeless and 
therefore changes over a length that is at least as great as that which char- 
acterizes the variation of the potential. For a particle moving in a potential, 
such a wavefunction describes a bound state. In a bosonic field theory it is 
describes the ground state. 

The anti-WKB approximation is not as straightforward as the WKB ap- 
proximation because it is not generally true that IV 2 ^! >> \VS\ 2 for bound 
state wavefunctions. To see this explicitly, assume the potential vanishes 
when r > a. Then for r > a, 



S = S,-Kr- I ^\n{r/r )- K = J-^, (1.3) 

where motion is in D spatial dimensions. Comparing the gradient and Lapla- 
cian applied to this expression, we find that the anti-WKB approximation is 
justified only when D = 1 and K = (or at least K ~ 0.) This is the case 
of a weakly bound particle moving in one dimension. (Recall that in one di- 
mension there is always one bound state in an attractive potential, no matter 
how weak the potential.) This version of the anti-WKB approximation has 
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been developed. [E] One obtains an expression for the energy of the bound 
state as a series of integrals over powers of the potential. The relative size of 
the n-th term is proportional to the n-th power of the parameter V^ma 2 /ft 2 , 
where Vq is the strength of the potential and a its range. We see here an 
expected contrast with the WKB approximation: a series of decreasing pow- 
ers of Plack's constant. Despite this, if the potential is such as to make the 
parameter small, the approximation succeeds. 

The straightforward case just described is quite different from the anti- 
WKB approximation considered in this paper, where we treat the general 
case D ^ 1, K ^ 0. We evade the "no-go" conclusion, above, by supplying 
an initial guess for S. Thus we write 

S(r) =F{t)+T{t), (1.4) 



where the variational seed F has the asymptotic form [L^, and is chosen to be 
an initial guess for S. The correction T is now determined by Schrodinger's 
equation: 

V 2 X„ + 2VF-V7; = ^ [V(r) - E n \ - V 2 F - |VF| 2 - \VT n ^\ 2 ; 

n 

r_i = o. (i.5) 

The index n = 0, 1 . . . labels the successive approximations to T and hence S. 
The initial approximation ignores the gradient of T; the next approximation 
uses the initial approximation for the gradient, and so forth. Of course, 
ignoring the gradient of T is different than ignoring the gradient of S, much 
of which appears in Eq. |1.5| in the term | V-F| 2 ; nevertheless, we still call this 



the anti-WKB approximation. Note that for each n we obtain a different 
approximation for the bound state energy, E n . 

In Section 2 we demonstrate that this sequence of approximations is for- 
mally convergent provided F is well chosen. It is useful to state what we 
find. Assume that our variational seed differs from S by a function scaled by 
a small pararmeter e: 

F(v) = S(Y)-eS l (v). (1.6) 
It follows that the exact T is eS\. In Section 2 we show that 

E n = E + e n+2 A n , T n = eS l + e n+2 R n . (1.7) 
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Here E is the true bound state energy, and the factors A n and R n are finite at 
e = 0. The convergence is formal because we have not given estimates for the 
n-dependence of A n and R n . If these factors increase with n faster than an 
exponential, the sequence does not not converge, but is instead asymptotic. 

The convergence estimates of Eq. |1.7| depend on the fact that the term 
T n _! appears quadratically on the right hand side of Eq. |1.5| . When F differs 
little from S, treating the gradient as a correction is justified. We again 
emphasize the difference between what is done in Ref. and what we do 
here. This is not a weak potential approximation, nor is it an expansion in 
inverse powers of h. We rely on the availability of a reasonable guess for the 
variational seed F. The guess must also be such that the integrals and Green 
functions we shall encounter can be computed. 

The wavefunction arising from our construction is nodeless. The reason 
is that if ip vanishes on the surface /(r) = 0, there will be a term ln/(r) 
in S. In the neighborhood of the node, this logarithm dominates S, and we 
must explicitly incorporate it into F if we are obtain a wave function with 
a node. Note that we must specify the surface / = 0, which is known only 
under special circumstances. We do not consider such cases here; our F's will 
be smooth, and we therefore limit ourselves to nodeless wavef unctions. This 
means that the energies E n are approximations to the ground state energy 
E. 

Like the WKB approximation, the anti-WKB approximation is nonper- 
turbative; it does not require the presence of a small parameter in the Hamil- 
tonian. Its major limitation is that it is restricted to the ground state wave 
function. But the anti-WKB approximation has this important advantage: 
Equation [L5] for T n is linear. We will see that it can be solved readily in many 
cases of interest, including particles moving in three dimensions in asymmet- 
ric potentials, and many-body problems like bosonic lattice field theory. In 
the latter problem the method can to extended to study some of the vacuum 
state correlation functions that are of central importance in field theory. All 
these possibilites are closed to the WKB approximation, which is generally 
unmanageable except for a particle whose motion effectively reduces to one 
dimension. It should be noted that the anti-WKB approximation is less 
useful for problems involving identical fermions because the Pauli principle 
makes the nodeless state unphysical for systems of more than two spin 1/2 
particles. 

In Section 2 we present a theoretical development of the anti-WKB ap- 
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proximation. We include the solution of the dynamical equation |1.5| , the 
convergence of the sequence of approximations, and several other matters. 

In following sections the anti-WKB approximation is applied to a number 
of problems to show that it works in increasingly complex situations. In 
Section 3 we study a spherically symmetric square well whose ground state 
is known by elementary methods. We find that E Q gives only 62% of the 
correct binding energy, but the next approximation, E\, gives 96%. We 
develop Green functions required for the application of Eq. [L5] to potential 
problems where the potential is nonspherical. In Section 4 we study long 
range potentials, with the hydrogen atom as a particularly simple subcase. 
We finish by applying the anti-WKB approximation to many-body problems: 
4 field theory in Section 5, and the helium atom in Section 6. Conclusions 
are presented in Section 7. 



2 Anti-WKB Equations 



The energies in Eq. |1.5| are determined by a general requirement. Consider 
the surface integral 



dA- 



'e 2F VT n _ 


= f dVV- 


e 2F VT n 




Jv 





(2.1) 



As S is expanded to infinity, the surface integral decreases to zero because of 
the exponential fall of the factor e 2F . Using Eq. |1.5| we obtain the eigenvalue 
equation determining E n : 







/ dVe 2F {|J [V(r) - E n ) - V 2 F - \VF\ 2 - |VT n 



-l 



(2.2) 



We can now assemble the equations to derive the convergence results 
|1.7| . The true ground state energy E is determined by Eq. |2.2| with the 
replacements E n — > E, T n _i — ► T. Subtracting equations, 







dVe 2F 
E„ = 



2m 

w 

E + 



(E„-E) + |VT n _ 1 | 2 -|VT| 2 

h 2 J dVe 2F [\VT\ 2 - [vr^xi 
2m 



(2.3) 



JdVe 2F 
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We now express the n-th order factors in terms of their limiting values and 
a deviation: 

E n = E + an, (2.4) 
T n = -eSi + r n . 

Then Eq. |2.3| becomes 

2 h 2 JdVe' F \VS 1 \' 
a ° = £ 2^ JdVe^ ' (2 - 5) 
ft 2 JWe 2F [2eV^ r Vr n _! - |Vr n _!| 2 ] 



2m / dVe 2F 



(n > 1). 

We next use Eq. [L5] and the analogous equation for T = eS\. Subtracting 
the equations 

V 2 r + 2VF-Vr = --^ao + e 2 | V^] 2 , (2.6) 

n 



2 777, 

V 2 r n + 2VF-Vr n = ~JT a n + 2eVS , r Vr n _ 1 - | Vr n _x| 2 . (n > 1). 



It follows that Eq. [O] holds for n = 0. Furthermore, if |1.7| holds for n, then 
from Eqs. |2~R ETH we see that it holds for n+ 1, and the result is established 



by induction. 

Energy £7 n requires T n _i for its computation. This information can be 
used to find a better estimate of the state energy from the expectation 

< ip\H\ip > 

n = ' ^ = 6XP( + n - lY ( 7) 

The error in ip is 0(e n+1 ); variational argument || gives 

£ n -E = 0(e^); ^f- |r = 0(e"). (2.8) 

Eq. |1.5| may be solved using the Green function satisfying 

LG(r,r') = 5(r-r'), (2.9) 
L = V 2 + 2(VF)-V. 
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Then 



T n (r) = JdV'G(v,r')S^[V-E n 



V 2 F- IVFI 



(2.10) 

It is illuminating to expand the Green function in terms of the orthonor- 
mal eigenf unctions of L, because the argument leads to the eigenvalue equa- 
tion |2.2j in a different way. 



dVe 



2F , 



■>kWk2 
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kl,k2- 



Then 



G(r,r') = E-^( r )^( r, )e 2F(r,) . 

k 



(2.11) 



(2.12) 



But there is a problem: L annihilates a constant function, so there is a 
normalized eigenfunction of L with eigenvalue zero: 



Ur) =N, N 



dVe 



2F 



-1/2 



Ho = 0. 



(2.13) 



The contribution of this zero mode to G is infinite. Nonetheless, T n in 
Eq. |2.10| is finite if the projection of the driving term onto the zero mode 
vanishes. The condition for that is the eigenvalue equation |2T2 . 

The eigenvalue equation for E\ can be simplified. We need the integral 



- / dVe 2F \T { 



e 2F T VTo\ + / dVe 2f T V- [e 2F VT l (2.14) 



dVV- 



The first integral on the right vanishes by the divergence theorem, and the 
second may be transformed using Eq. |1.5|. 



- J dVe 2F \T \ 2 = J dVe 2F T D = J dWe 2f(r) D(r)D(r')G(r,r'), (2.15) 

(2.16) 



where 



D 



V 2 F 



IVFI 2 . 



The quantization condition for E\ reads 

= J dVe 2F D + ^Eo-E^ + J dVdV'e 2F ^ D(r)D(r r )G(r, r'); 
El = Ev I H * JWe 2F( ^(W)g(^) i (217) 



2m 



J dVe 2F 
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When V(r) depends on a single coordinate, all differential equations in- 
volved in constructing the Green function can be solved. Consider the case 
where V, F and T depend only on the radial coordinate, as happens in Sec- 
tion 3. Then the Green function satisfies the ordinary differential equation 

Lg(r,r') = -^5(r-r>), (2.18) 

L- — + 2(- + — \ — 

dr 2 \r dr J dr 

Rather than expand in eigenfunctions, we construct g from the two solutions 
of LR = 0, which are known. 

r dr'e' 2F ^ 

R + (r) = 1; R-(r) = / . (2.19) 

Jro r 

The Green function is 

g (r,r') = R + ( r> )R-( r< )/r' 2 W(r') = -R + ( r> )R'( r< )e 2F ^ , (2.20) 
where W is the Wronskian 

dR+ dR~ e - 2F ( r ) 

W(R\ R-) = ^-R- - R +< ^- = (2.21) 
dr dr r l 

Now 

roo [ O777 

T n (r) = j o dr'r' 2 g(r, r>) | ^ [V - E n \ (2.22) 
-F"-^-F'-{F'f-{T' n A . 

' / r' 

The eigenvalue equation emerges here, not through a zero mode, but through 
the bad asymptotic behavior of R~. Near r = 0, R~(r) ~ e~ 2F ^ /r, and 

Ur) ~ ^— jf drW F ^^ [V - E n ] (2.23) 

-F"-lF'-{F'f-{T' n _A . 

To keep T n finite at r = 0, the eigenvalue integral must vanish; when it does, 
T n (0) = 0. 
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3 Spherical Square Well Potential 

We begin our study of the anti-WKB approximation by applying it to the 
case of the spherically symmetric potential 



V(r) 



-V Q (r < a) 
(r > a) 



(3.1) 



The radial Schrodinger's equation for s-wave bound states can be solved, and 
the energies are determined by the equation pfl 



ka cot ka 



-Ka: 



ka 



! 2ma 2 



V - (Ka 



2. 



E 



h 2 K 2 



There is no bound state unless Vq exceeds the threshold strength 



(3.2) 



esh 



71 



Ci; 



d = 1.0, 



(3.3) 



2ma 2 \2 

and for Vq just above threshold strength, the energy of the ground state is 



E 



ma 



{Vq — [V }thresh) 2 C2', 



Co = 1.0. 



(3.4) 



To compute E we use the eigenvalue integral |2.23| . To get some notion 
of the dependence of E on the variational seed, we try three of them, each 
of which depends on a single parameter tq\ 



FJr) 



Fi(r) = —Kr — ln(r/ro + 1) 
F 2 (r) = -Kr - In \/(r/r ) 2 + 1 
fl-(K + ±)r, (r<r ) 



(3.5) 



-Kr — ln(r/ro), (r > ro). 



(F3 is chosen so that the function and its derivative are continuous at r = ro.) 
Using these variational seeds, we obtain the results summarized in Table 2.1. 
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F(r) 


r /a 




c 2 


S(r) (exact) 




1.0 


1.0 


F 1 (r) 


0.638 


1.14 


0.127 


F 2 (r) 


0.923 


1.09 


1.37 


F 3 (r) 


1.38 


1.03 


0.826 



Table 2.1. The optimum values of r$/a and the coefficients C\ and C 2 for 
the three choices for F(r) . 

Note that Eq = £ is a variational energy; in each case the parameter r /a 
has been chosen to minimize E . We see that the results are sensitive to F. 
F 3 gives a threshold potential strength that is about 3% high, and it also 
gives the best estimate for the bound state energy near threshold. We use 
F 3 in the rest of Section 3. 

To explore further we choose a potential strength well above threshold: 
2ma 2 Vo/fr 2 = 3.0 The bound state energy is given by Eq. |3.2| to be E = 
— 0.0613k 2 /2ma 2 . We find that at this potential strength E is minimized 
for r /a = 1.505 (in the threshold calculation the value was 1.37). We find 
E = -0.0379h 2 /2ma 2 . A measure of the (mediocre) quality of this result is 
the ratio E /E = 0.619. 

We next compute E\. For this one degree-of- freedom problem, we use 
Eq. |2.2| because a simple formula for dT n j dr is available: 



dT n e~ 2F ^ r , 2 f 2m 



dr r 2 

2 



F"{x) - -F\x) - [F\x)} 2 - [T'^x)} 2 . (3.6) 



x 



We now obtain much improved results: E% = — 0.0589fo 2 /2ma 2 , and E\/E 
0.960. Almost 90% of the error in Eq has been removed in E]_. 
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One of the properties of the anti-WKB approximation we have empha- 
sized is that it can be applied to potential problems where V(r) has a general 
dependence on r. This statement is qualified by the requirement that F(r) 
be such that Eq. [2.9| can be solved for G. One case where this can be done is 
when V is nearly spherically symmetric. Then we choose F to depend on r, 
relying on the T n to supply the nonspherical corrections. Note that the anti- 
WKB method is more flexibile than the method of separation of coordinates 
in that it is not necessary for surfaces of constant potential to exactly fall on 
surfaces of constant coordinate. Approximate coincidence suffices. Below we 
use separation of coordinates to construct G in spherical coordinates. 

When F depends only on r, we write G in the form 

OO 7YI — £ 

G(r, r') = $>(r, r') E 0)^(0', </>'), (3.7) 

£=0 m=-£ 

where the radial Green function satisfies 

d 2 ge(r,r') , (\ , dF\ dg e (r,r') £(£ + !). 1 

)( 1 g e (r 1 r) = —5{r-r). (3.8) 




This Green function is 



9i\r,r 



RtirJRtirJtPW/Ct, (3.9) 



where R + ' satisfy the homogeneous equation 



^ +2 (l + £)f£_4£ + i)* = . (3.10) 
dr z \r dr J dr r z 



The factor Ct comes from the Wronskian W{R£,Rj) = C e e- 2F ^/r 2 . 

For £ = 0, Eq. 3.10| has the same solutions we found for the case of 



spherical symmetry, and we can choose R ' = R + ' of Eq. |2.19| . We again 



encounter the eigenvalue integral 2.23 with this change: Only the s-wave 



projection of V enters into the determination of Eq. This is in keeping with 
our assumption that V is nearly spherically symmetric. The higher E n 's 
receive contributions from all of the angular projections of V. 

The solutions of Eq. |3.10| for £ > require further analysis. We first note 



that near r = they behave like / orr x , and we must choose R t to have 
the nonsingular behavior r . At large r, solutions of Eq. |3.10| behave like 
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those of Eq. [2.19| because the last term in Eq. |3.10| is unimportant at large 



r. We must choose the solution behaving like a constant so that R + remains 
finite at large r. 

In the case of F 3 (r), Eq. p.lQ has solutions that are familiar functions. In 
the interval r < r , these are confluent hypergeometic functions: 

R a = /$(£, 2£ + 2, 2Kr + 2r/r ), (3.11) 
R i2 = r -^- 1 $(_^ _ l ; -21, 2Kr + 2r/r ). 

In the interval r > r$ they are modified spherical Bessel functions: 

Res = V^e Kr I e+1/2 (Kr), R u = V^e Kr K e+l/2 (Kr). (3.12) 

We then have the solutions 

1 \ aR e3 + /3R U , r>r ; e \ R e4} r > r . 1 ' ' 

The coefficients are fixed by demanding that the functions and their first 
derivatives be continuous at r = r$\ they may be expressed in terms of 
Wronskians evaluated at r = r . 

When a particle moves under the influence of a highly nonspherical po- 
tential, the Green function should be constructed in coordinates chosen so 
it is reasonable for F to depend only on one of them, £. As an example, 
consider a potential well in the shape of a right circular cylinder of diameter 
D and length L: 

y( T ) = t-V \z\ < L/2, Vx 2 + y 2 < D/2; (3J4) 
1 otherwise. 

If L ~ D, this is an example of an approximately spherical potential, and it 
is appropriate to use the Green function constructed above. Here we consider 
the case of a long thin rod: L » D. Now it is manifestly inadequate to 
take the surfaces of constant F to be spherical. These surfaces ought to be 
figures of rotation about the z-axis in the form of elongated cigars: prolate 
spheroids. This suggests use of the prolate spheroidal coordinate system: 



11 



1 sin 9 cos ( 



r 



1 < £ < oo, 
dV = 



£ 2 — 1 sin 9 sin 0, 

z = — E cos 9; 
< (A < 2tt; 



(£ 2 -cos 2 #) 8va.6d£d6d4>. 



(3.15) 



< < 71, 

2 



Surfaces of constant £ are ellipsoids of revolution about the z-axis, the ellipses 
having foci at z = ±r /2. If r is appropriately chosen (by a variational 
minimization of E ), the shape of the ellipsoids can be made to resemble 
that of the rod for £ ~ 1; for large £ the ellipsoids approach spheres of radius 
ro£/2. An obvious choice for F is 



F(0 = -tfr £/2-ln£. 
In prolate spheroidal coordinates, Eq. |2.9| takes the form 



d_ 



(^ 2 -D- 

dE, 

+2(e - 1) 



+ 



i d 

sin 6> (96 1 



Sm9 d9 



+ 



1 1 

+ 



£ 2 -l sin 2 fl 



<9 2 G 



dF dG 



The solution has the form 



r 



(3.16) 



-5(C - O s ( cosd ~ cos #')£(</> - </>')• (3.17) 



<=0m=-< 



where the radial Green function satisfies 



(^-l)?§P + 2 



,9 N rfF 



dg 



trn 



The radial Green function is 



)Y, 



+ 1) + 



m 



(3.18) 



9 in 



— ^-O- (3-19) 



(3.20) 
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Cim appears in the Wronskian W(E + , H~) = CV m e~ 2F /(£ 2 — 1). The func- 
tions Ef m satisfy the homogeneous equation obtained by removing the delta 
function source in Eq. |3. 19 . 

The general method of constructing the Green function appropriate to 
the rod-shaped potential well is now clear. We discontinue discussion of the 
functions Ef m except to point out that the differential equation they satisfy 
has three regular singular points and an irregular singular point at infinity. 



4 Long Range Potentials 

When a potential is long range and behaves at large distance like r~ p , 1 < 
p < 2, there is no longer a threshold strength for the bound state. The 
anti-WKB approximation reproduces this result, and it does so in lowest 
order, when n — 0. We demonstrate this by showing that with the choices of 
F given by Eq. [3.5| we can adjust tq so that Eq=0 no matter how weak the 
strength of the potential tail. A futher adjustment of r$ will then produce a 
negative (bound state) energy. 

We consider the case of F 2 . For the potential V /r p , (r > a), the eigen- 
value equation implied by Eq. [2.23] takes the form, when n = E Q = 0: 



2mV f°° r 2 p dr 2 f°° r 2 dr 



°" h 2 la r 2 + r 2+3r °Jo (r2 + r )3- ^ 



This may be written 



mV rl- p 2vr oo/oya+3-p (_i)fc 3?r 



hf 2sin7r(p- l)/2 ^ \r J fc + (3-p)/2 16r 

For 1 < p < 2, there is some large value of ro for which this equation is 
satisfied, no matter how small Vq may be. There must be a nodeless bound 
state. 

When p = 1, we come to the classic case of the Coulomb potential. 
At this point, the asymptotic behavior of the wavefunction changes, and in 
three dimensions the coefficient of the ln(r/r ) term depends on the potential 
strength and bound state energy. Since this coefficient is no longer known, 
the obvious response is to treat this coefficient as an additional parameter, 
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T\. Thus, for example, F 3 is generalized to 



F,(r) = l ri - {K+ n^ ( r<r °) ; (4.3) 
t — Kr — T\ ln(r/r ), (r > r ). 

But now the determination of the variational parameters is obvious, because 
with r\ = 0, F3 = — Kr = S. Therefore T{r) = 0, and we immediately have 
the hydrogen wavefunction, with E = E. The same result is achieved with 
Fi and F 2 . 

It is fortuitous that the hydrogen wavefunction occurs among the nat- 
ural choices for F. Still, it is nice that the Coulomb potential is so easily 
encompassed by the anti-WKB approximation. 



5 4 Field Theory 

The anti-WKB approximation can be applied to problems having many de- 
grees of freedom. We illustrate this using the case of 4 field theory in one 
spatial dimension. In the continuum the Hamiltonian is 



H 



l r 1 1 

dx \ 2 CV(3;) + 2 



d(f)(x) 
dx 



1 c 2 

+ -j^{m 2 + 5m 2 )4> 2 {x) 
+ 24^4 



(5.1) 



We identify fields at x = and x = L. To treat the system by the anti-WKB 
method we must deal with a discrete set of degrees of freedom. We therefore 
divide the line into N segments of length a; N = L/a. Each segment is 
represented by a dimensionless lattice coordinate <p k , k = 1, . . . , N. The 
lattice Hamiltonian is 



H 



he" 



1 d 2 



2 901 



+ ^(0 fc+1 - fc ) 2 + \ (y) (™ 2 + bm 2 )<p 2 k 



Xhca 
24 



4 



(5.2) 



Putting the field theory on the lattice introduces a short distance cutoff 
a, or equivalently a large momentum cutoff h/a. This removes the notorious 
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divergences of continuum quantum field theory. However, the divergences 
still lurk and reveal themselves when we approach the continuum limit by 
taking a — > 0, iV — > oo, with Na = L fixed. We want to be able to take the 
limit, of course, since for us the lattice is only a computational device. In our 
case we find that as we take a — > we approach a continuum theory having 
infinite mass. (We show this below.) This pathology disappears when we 
include an appropriate mass counterterm in the Hamiltonian: 




This counterterm is the lattice version of the very counterterm that must be 
included in the continuum Hamiltonian. At small a, A grows logaritmically 
to keep the effective mass finite. 

l h 2 , . 

A ~ - In 5.4 

7r m z c z a z 

In one spatial dimension 4 field theory is superrenormalizable. For 
this system, the explicit counterterm in |5.3| suffices to remove all a —>■ 
divergences in "physical" entities like the correlator 

1 N 

< O|0(O)0(x)|O >~ he < 0|- <Pk(p k+x /a\0 > . (5.5) 

fc=i 

Unfortunately, the groundstate energy E that has figured prominently in the 
anti-WKB approximation is not "physical" in the sense used here. Even free 
field theory, A = 0, includes zero point energies for each degree of freedom, 
and when summed these diverge as a — > 0. We will cope with E as we go 
along. This point aside, Eqs. |5.2| , |5.3| define a lattice representation of con- 
tinuum </> 4 field theory to which we can apply the anti-WKB approximation. 

A central issue for this many degree of freedom problem is finding an F 
that is appropriate, and for which the multidimensional integrals and Green 
function can be found. There is just one choice for which this is easy: 

F = ~\ E ^WfcA 2 = -~<W, (5.6) 

k\k2 
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with M a real symmetric matrix. Consider the eigenvalue equation for Eq. 
= J (Ih#)e-^{]T[(0 fc+1 - fc ) 2 + (ca/h) 2 (m 2 + 6m 2 )tf k 

k 

+ (\hca 2 )<f)i/ 12] - 2E a/hc + Tr(M) - (f)M 2 (j)}. (5.7) 

Every term in the brace involves powers of so the integral can be evaluated 
by taking derivatives with respect of a k of the generating function 



I\[a) = J (Hdcj)) exp(— 0M0 + iac, 



^a k (f> k . (5.8) 

k 



This generator, in turn, can be evaluated easily by a change of coordinates. 
Let the eigenvectors and eigenvalues of M be v q and \i q . 



Mv q = fi q v q , 
The new coordinates are 

k 



(q = l,...N). 



(5.9) 



(5.10) 



In these coordinates I\ becomes a product of N independent Gaussian inte- 
grals, and 

' • 1 4 (5.11) 



. exp [ — clM 1 a ) . 
VdrtM V 4 J 



Eq. |5J becomes 

= h(0){Y,[Mkk - M^+i + (ca/h) 2 (m 2 + 5m 2 )M kk 1 /2 
k 

+ (Xfica 2 )(M k - k 1 ) 2 /lQ\ - 2E a/hc + Tr(M) - Tr(M 2 M~ 1 )}. (5.12) 

M fclfc2 must be a function of fci — k 2 of period iV owing to the homogeneity 
and periodicity of the Hamiltonian. It can be written as a discrete Fourier 
transform. 



1 



N 



M^ k2 = -Y,^\p)^V 
JV P =i 



2Tiip{k\ — k 2 ) 
N 



(5.13) 
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With this the eigenvalue equation decouples into modes. 



I p 



u{p) 



cos ■ 



2np 

w 



(ca/h) 2 (m 2 + 5m 2 



1 E 



N V w(p) 



Tic 



(5.14) 



+iV(Aftca 2 ) 

The cu(p) are parameters in M that are fixed by minimizing Eq. We find 

u(p) = 



ca\ 2 , _ . os A/ica 2 . . 9 
(m 2 + 5m 2 ) + H — A + sm 2 



n J 



up 

~~N 
1 



1 E 



(5.15) 
(5.16) 



Eq. |5.16 determines A, and it is enlightening to examine its solution near 
the continuum limit when we drop the counterterm Sm 2 . Following Eq. |5.4| , 



A 



-In ; 

7r A/ica 2 A' 



A 



In 



4/4 
+ In In 



7T 



Xhca 



Xhca? 



(5.17) 



Therefore, the square of the effective mass in uj{p) becomes, in the limit of 
small a, 



m 2 + 



Aft 3 A 



\h A 



4c + < 5 ' 18 ' 

This shows that we approach a continuum theory having infinite mass. On 
the other hand, when we do include the counterterm, Eq. |5.16j has the solution 
A = A ; u)(p) = LOo(p). Then the mode function ui{p) never changes, no 
matter what the lattice spacing or coupling strength. 
The energy Eq is 



Eq 

T 



he 

2a? 



N 



P =i 



i ^ 
m c 

~h~ 



Xh 

m 2 c 



A 2 

0' 



(5.19) 



The sum on the right is the zero-point energy of the degrees of freedom. It 
is present in free field theory, and its contribution is quadratically divergent 
in a -1 . The last term, proportional to the dimensionless coupling Xh 3 /m 2 c, 
diverges like (lna -1 ) 2 . In perturbation theory it arises because anomalous 
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combinatorics spoil the cancellation of divergences in self-energy loops and 
A . This failure to cancel occurs only at order A of perturbation theory, so 
Eq already exhibits all the terms that diverge as a — > 0. 
To compute E\ we need the Green function satisfying 

[w ~ 20M ^) 00 = ? 6{<pk ~ ^' (5 ' 20) 

in matrix notation. We use Eq. |2.12| to construct G, which requires us to 
construct eigenfunctions of the operator on the left. We again use coordinates 
% then the equation for the eigenfunctions and eigenvalue is 



N 

E 

q=l 



d(r) 



d 2 d 
tL 2y?tf— 



q\2 



ip = flip. 



(5.21) 



This equation may be solved by separation of variables. 

N 
9=1 

The functions H s are Hermite polynomials satisfying 

d 2 Hjx) dH 8 {x) . 

— 2x ; h 2sH s (x) = 0, 



dx 2 



dx 



/oo 2 
dxe~ x H sl (x)H S2 (x) = 5 SuS2 ^n2 Sl (s 1 \). 
-co 



(5.22) 



(5.23) 



There are N indices on our wave function, and the eigenvalue is a linear 
function of them: 



N 



q=l 



The Green function is 
1 



G 



-2 ]e 



-4>'M<j>' 



' dz -J-J- I fig 

9=1 



X 



(5.24) 



(5.25) 
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The purpose of the integration over z is to produce the denominator lf/ik in 
Eq. 2.12| . The sum over s q is given by Mehler's formula. |7j] Reverting to the 
original variables, 



1 Jdet(M) 



dz 



- z^Jdet{l-z 2M ) 
, Mz 2M 



exp 



M 



1-z 

M 



2M 



y 2M 



+ 20 



, Mz 



z 



2M 



(5.26) 



As we expect from the discussion of Section 2, this Green function does 
not exist because \i = is an eigenvalue. In Eq. |5.26| the divergence appears 
at z — 0, and it has been regulated by a cutoff at z — e. However, when 
computing T n , the coefficient of the 1/z factor is 



-<j>M<t> 



D0; 



(5.27) 



D 



XTica 2 
12 



E 

. k 

2aE, 



y 2 + 



ca\- 



(m + 5m 



2\a2 
k 



Tie 



- + Tr(M) - <f)M 2 (f) - ( 



Setting this to zero, we obtain a finite T n as e — > 0. Thus we again en- 
counter the eigenvalue equation as a consistency condition for the solution 
of Schrodinger's equation. 

The dependencies on the fields in Eq. |5.26| are Gaussian, so we can use 
Ji to evaluate the generating function 



I 2 (a,/3) 



;JJ d(f)d(f)')G((j), 0') exp(-0M0 + ia<p + z/30') 

n N/2 



cxy> [-~aM~ l a- -(3M~ X (3 



exp 



~-a[z M M- l ]p 



'det(M) 

i\ r 1 dz 

2 

By differentiation of I<i we obtain all the contributions in 

hcf(Ild(t)d(t)')e-'f >M ' t> G(<j), <j/)D{<j>)D{^ 



(5.28) 



Ei — E 



2a 



SW)e- 



(5.29) 
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The result is 



where 



A, 



E!-Eq m 2 c 3 f\h 3 \ A 2 
L " h~ V^Fc) 384 : 

mca\ 2 f 1 dz 1 



Jo z N £-t V Jkl ' k2 



Using Eq. |5.13 



I N z ^o(p) 



(z M M~ l ) k k — — ~T~S exp 
n ^ w (p) 



2irip(ki — A^) 
TV 



(5.30) 



(5.31) 



(5.32) 



fmca\ 2 ( 1 \ 3 ^ 



1 



X 



Pi,P2,P3 ^(PlVofeVofeVo^l +P2 +P3) 
1 

[W (pi) + Wofe) + ^o(P3) + Wb(pi +P2+ Ps)] ' 



(5.33) 



It is straightforward to verify that A 2 remains finite as a —>■ 0. 

What emerges after all our work appears to be perturbation theory in 
"old-fashioned" form. This is easily demonstrated. We have E\ — Eq = 
0(A 2 ) = O[(dT /d(f)) 2 }. Hence, T = 0(A). Then use 



T n -T n _ 1 = J{m<p')G{<p,<p') 



2a 

he 



(E 



n-l 



'<9Tn-i* 



(5.34) 



E n — E. 



S(Ild<P)e 2F [(dT n _ 2 /d<P) 2 - (dT^/dePY 



n-l 



SW)e 



2F 



to prove by induction T n — T n _i = 0(A n+1 ) = E n — E n -\. This result is not 
surprising. Our F is the S of free field theory, and we compute corrections 
to that. This is the program of perturbations theory. 

No one wants to use the anti-WKB approximation to develop perturba- 
tion theory. What we learn from our exercise is that the calculations can be 



20 



carried out, and they lead to a sensible and familiar result. But it is possi- 
ble, within the anti-WKB approximation, to contemplate alternatives that 
transcend perturbation theory. An obvious generalization of F is 



N r 1 1 

2 W) 2 - f q (v q T 



E 

q=l 



(5.35) 



We use plane wave coordinates so that all degrees of freedom are coupled. The 
mode eigenvalue equation is unfamiliar, and Mehler's formula is not available, 
so we must work harder. A computer would be required to assemble G. 

A payoff is possible when a seed like F is used. Recall that at \h 3 /m 2 c ~ 
10, 4> 4 field theory makes a transition to a state of broken <p <-> — symmetry 
in which the field has a vacuum expectation value: < O|0(x)|O >= O 0- @ 
From this fact alone we see that F must be inadequate at strong coupling 
because with F, 0o — 0. But in F, as the variational parameters jl q vary 
smoothly with A, a vacuum expectation value can develop beyond a critical 
coupling. 

A trick attibuted to Feynman applies to this problem. Extend H by 
adding a term £(l/iV) Y,k ^k^k+d- The extended Hamiltonian remains ho- 
mogeneous and periodic, so the techniques we have developed continue to 
apply. E depends on £, of course, and by first order perturbation theory (in 
£), dE/d£(£ = 0) is just the correlator of Eq. |5J| In this formula we would 
substitute E n for E. 



6 Helium 

Because of the Pauli principle, helium is the only multi-electron atom in 
which a symmetric spatial wave function is physical. The nodeless state we 
construct is the groundstate of helium. 

The Hamiltonian we study includes only the potential terms represent- 
ing the electric forces acting on electrons moving around an immobile nu- 
cleus. The ground state energy of the simplified Hamiltonian therefore differs 
slightly from that of physical helium where magnetic, recoil and relativistic 
corrections are present. The Hamiltonian is 

z r 1 r 2 ri — r 2 
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We use dimensionless coordinates, so the eigenvalues of H must be multiplied 
by me 4 /ft 2 . 

A standard textbook variational computation with this Hamiltonian as- 
sumes that the wavefunction is a product of exponential (hydrogenic) fac- 
tors. || The resulting energy < H >= — (27/16) 2 = —2.848, is very close to 
the experimental groundstate energy, -2.90, and we adopt it as a standard for 
the problem. Since we have the the tools to deal with many-body problems 
when F is a quadratic form, we use that: 



Our choice of F corresponds to a Gaussian first guess for the helium wave- 
function. We compute Eq, choosing Kq and K\ to minimize it and find 



We obtain only 0.816 of the "standard" binding, which means that the anti- 
WKB correction is substantial. 

An interesting feature of our result is that Ki is so small. The repulsion 
between electrons enhances the probability of finding the electrons on oppo- 
site sides of the nucleus, but not by much. In fact, if we set K\ = 0, the 
binding is decreased by less than 1%, and Eq = —2.301. We use this version 
of F when calculating E% because of the resulting simplifications. The com- 
putation of E\ now parallels that of Section 5. The main new feature is that 
the terms in the potential are Coulomb, not polynomial. When we use gen- 
erating function I 2 , we integrate over parameters to generate the Coulomb 
terms, using the relation 



F = - l -K Q {r 2 1 +rl)-K l v r v 2 . 



(6.2) 



E = -2.324; 



K 



1.549 



K x = (0.09238)A^ ( 



(6.3) 




(6.4) 



Summarizing our results: 



E = -2.301; Ex 
E 



0.808 



2.707. 



< H > 
Ex 



2.848; 



0.951. 



(6.5) 



< H > 



<H> 
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7 Conclusions 



We have proposed an approximation for quantum mechanics and field theory 
that assumes the relative importance of gradient and Laplacian terms is 
the reverse of that in the WKB approximation. As with any such scheme, 
usefulness is an important issue. We have explored usefulness by applying the 
anti-WKB approximation to several problems. The examples of the spherical 
square well and the helium atom provide interesting numerical results. In 
these cases we chose a simple variational seed F and found that the initial 
(variational) estimates of the ground state energy were incorrect by several 
tens of percent. The first nontrivial anti-WKB correction removed three 
quarters of the error, or more. 

In practice, the WKB method is difficult to apply to problems with many 
degrees of freedom. The anti-WKB method does suffer from this fault. With 
a quadratic choice for F we developed the formulas necessary for the treat- 
ment of the rather different problems of 4 field theory and the helium atom. 
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